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I. INTRODUCTION 

B: 

Colloidal suspensions contain mesoscopically large particles from the nm to /im size regime, the colloids, together 
with small solvent and solute molecules. These molecules are orders of magnitude smaller than the colloids, and 
still, they heavily influence the interactions between the colloidal particles. Such microscopic species also significantly 
i^i - affect the thermodynamics of the system. However, they often cannot be directly visualized and most experimental 
techniques such as small angle X-ray or neutron scattering, probe the colloidal degrees of freedom only. This stems 
t-H , from the wide separation of characteristic time and length scales between microscopic and mesoscopic degrees of 
freedom, which therefore suggests to develop effective approaches for the colloidal particles, integrating over their 
, microscopic counterparts, that follow adiabatically. 

The focus of this manuscript will be on the behaviour of charge-stabilized colloidal suspensions Q, 0, In the 
\ subsequent analysis, we will further restrict ourselves to a mean field treatment where the correlations between 

■ microions are discarded. The resulting Poisson-Boltzmann (PB) description is valid in aqueous solvents under usual 
\Q ' conditions of temperature and for monovalent microions, because existing colloids do not allow to reach the high 

Coulombic couplings that are required to observe deviations from mean- field 3] . In the following, we therefore focus 
on the more special case of charged colloids in a simple 1:1 electrolyte, adopt the PB approach, and treat the micro-ions 
as point-like particles dissolved in a continuum of solvent molecules. 

Even within such a simplified framework, the solution corresponding to N c interacting colloids is a difficult problem 
from a numerical point of view 0,0- In principle, one has to compute the density distribution of microions in space 
Pmicro(^) for any given colloidal configuration. First one has to solve the PB equation around a large collection of 
iV c colloids. One can then compute the corresponding stress tensor to obtain the force felt by each colloid, from 
which the colloidal configuration of the next time step can be found. This defines the loop that has to be iterated 
- a numerically rather demanding task 0, 0, ^| . It is therefore still of interest and common practice to map the 
multi-component Hamiltonian of the charge-stabilized colloidal suspension onto a one-component model (OCM). In 
the OCM the colloidal particles interact via effective pair-potentials u e g 0, - effective in the sense that these 
pair-potentials reflect not only the true pair-potential between the "naked" colloids, but also the indirect effect that 
all the micro-ions and solvent molecules have onto the interacting colloids. As an advantage of the OCM, the number 
of degrees of freedom is significantly reduced: If each colloid has Z counter-ions, then a system with N c colloidal 
particles has at least ZxN c interacting species (not to speak of any additional salt ions or the solvent molecules), 
which in the OCM are reduced to just N c interacting colloids. 

The effective pair-potentials between charged colloids are usually Yukawa like (screened Coulomb, see Eq. (|TJ 
below), but with effective charge and screening parameters whose dependency on the colloid density makes the whole 
pair-potential density-dependent. There is no rigorous way how to "derive" these effective parameters from first 
principles. However, a large number of rather sophisticated recipes to determine effective parameters can be found in 
literature but no matter how sophisticated a scheme is, an heuristic element within these theories can never be 

avoided. In section [n] we will briefly describe two rather simple approaches to compute effective Yukawa- parameters, 
in the following referred to as the "PB cell model" (see e.g. ](|) and the "Jellium model" 0. 
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Combining various numerical simulation techniques (Monte Carlo, integral equation, and A^-body mean-field com- 
putations), we here investigate the performance of these effective Yukawa potentials and the two effective-charge 
theories, the PB cell model and the Jellium model. We will perform these tests in three different ways: (i) we first 
investigate if the form of the effective Yukawa potential is able to reproduce the correct distribution of forces felt by 
the colloids in various configurations. Such a comparison, that will be addressed in section lTTTl requires the solution of 
PB problem around N c colloids. One may anticipate that in situations of high salt content, the Yukawa approximation 
will be operational, while the opposite low salt regime deserves more attention, (ii) We then study the compressibility 
X of the suspension as a function of the colloid density. Both the PB-cell and the Jellium model do not only provide 
us with effective parameters for the Yukawa pair-potential, but directly predict the pressure and compressibility of the 
suspension. The Kirkwood-Buff relation now allows to check the thermodynamic consistency of the effective-charge 
models; it relates x to S(0), the infinite wavelength limit of the colloid-colloid structure factor. While \ follows 
directly from the two effective charge models, 5(0) is related to the structure of the suspension which can be obtained 
from the colloid-colloid pair-correlations calculated using the effective Yukawa forces of the OCM. As these Yukawa 
forces again require the effective parameters of the two effective charge models, the Kirkwood Buff relation can be 
used for a stringent test of the consistency of the effective charge models under scrutiny. Such a consistency check 
is performed in section llVl (iii) The most direct way to check the quality of an approximation is to compare it to 
the results of more rigorous approaches. In the third main section of this paper (section 01 , we compute osmotic 
pressure, compressibility and pair-correlations of the suspension for the two effective-charge models and compare the 
results to primitive model calculations. 



II. TWO MODELS FOR EFFECTIVE CHARGES AND SCREENING LENGTHS 

We start by outlining the two effective charge models that we here are concerned with. The idea behind using 
effective parameters to incorporate effects of non-linear screening in a pair-potential based otherwise on linear theory, 
is explained and reviewed in Ref. psj . 

Although the effective potential has an unambiguous definition, there is no rigorous operational route to construct 
this object and it is common belief, under scrutiny here, that a Yukawa form |l|, LJlil 

„ , x 7 2 x /exp(K off a)\ 2 exp(-K cff r) 

(3u cS (r) = Z cS X B — (1) 

V 1 + «cff a J r 

provides a reasonable description. Here a is the radius of the colloid (assumed to be a sphere), Xb is the Bjerrum 
length, kT — f5~ x is the thermal energy, while Z e g and n~g are the effective charge and screening length. Such a 
"DLVO"-like expression would accurately reproduce the large distance interaction of two colloids immersed in a 
salt sea HLE1E| but it should be kept in mind that it has to fail at short distances Q . One should also keep in mind 
that even if one had the perfect effective potential for two colloids, one can not be sure that it is also the appropriate 
effective pair-potential for a concentrated suspension of colloids, because in general these effective pair-potentials 
cannot be superposed. Many-body interactions between the colloids must be taken into account as they can have 
a considerable impact on the structure of the suspension |l9j . We also remark, that in the salt free case, the vali dity 
of is not obvious. Other limitations of the concept of effective pair-potentials are discussed elsewhere PII2,SEJ- 
In appendix ^] we recapitulate the essentials of the two empiric models (Jellium and PB-cell model) used here to 
determine Z e g and n~g in the Yukawa pair-potential. As emphasized in the introduction, effective parameters can 
also be used to predict the osmotic pressure of the suspension. We define the osmotic pressure from the pressure P 
and the reservoir pressure P re servoirj as 

II P -^reservoir P 2>C s kT , (2) 

where c s denotes the (monovalent) salt concentration in the reservoir. Furthermore, the osmotic isothermal compress- 
ibility x is defined via Eq. @ as 



* = — 



(3) 

T,salt 



where the derivative with respect to the colloid density p c is taken at constant salt chemical potential (i.e. constant 
c s ). Applying these two definitions, the osmotic pressure of the suspension can be related to the effective parameters 
(see appendix [XJ 

47rABa 2 /3n micro = n 2 s a 2 - K 2 a 2 (4) 
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with K 2 = 8tt\bc s . The suffix "micro" is explained further below. The compressibility then reads 



"off" 



3A_b drj 



(5) 



in which r\ — 4wp c a /3 is the colloid volume-fraction. In these expressions n^ s a is given either by the PB-cell 
expression, Eq. (|A3|) and l|A4() . or the corresponding expression for the Jellium model, Eq. (|A2|) . 

While the Jellium model is implicitly based on the assumption that the colloid-colloid pair distribution function 
g(r) ~ 1, the PB cell model takes the opposite view and assumes a g(r) corresponding to a crystalline colloidal 
configuration. This latter assumption rests on the observation that the repulsively interacting colloids arrange their 
positions such that each colloid has a region around it which is void from other colloids and which looks rather similar 
for different colloids. In other words, one model assumes a regular arrangement of the colloids in the suspension and 
thus a rather solid-like structure, while the other model assumes no structure at all - a situation typical of fluids at 
low density. In that sense, the two models chosen are complementary to each other. 

The effective charges as a function of the bare colloid charge i?bare go to Zb are for low bare charges, and to a 
saturated effective charge Z^ for high charges, and can therefore be roughly approximated by 

7 (7 \ ~ I ^barc ^bare ^ ^cff ( a\ 

^cffl^barcj ~ \ ysat 7 ^ 7S at V ) 

I ^cff ^barc £ -^ e ff 

Except in section we will mostly concentrate on the saturated regime that is most suited to describe colloidal 
suspensions [34| (for a discussion of the crossover behaviour between low and high bare charges see [371 for symmetric 
electrolytes and [38| for 1:2 and 2:1 electrolytes). By choosing in all our calculations a Zb are large enough to ensure 
saturation of the effective charges, we can thus get rid of the additional parameter Zbare- 

We have plotted the effective saturated charges as a function of rj in Fig. JSJ of the appendix. One can identify 
a regime in rj where salt-ions dominate the screening, leading to effective charges that are practically independent 
of the volume-fraction. At other values of rj however counter-ions outnumber the salt-ions; this is the regime where 
screening is dominated by the counter-ions. As the number of these counter-ions depends on the number of colloidal 
particles, the effective charges become ry-dependent when the counter-ions dominate the screening. One may estimate 
from Eq. (|A2f) of the appendix the volume-fraction rj* where one may expect to find the crossover between both 
screening regimes, 

* k2q3 fn s 

V* - Trn- • (7) 

We remark that this threshold value rj* was derived within the Jellium framework. We have no clear definition like 



this for the cell model and we have determined the rj* for the cell model empirically (see section IIVB(l . Empirically, 
as will be discussed in the section TlV Bl the crossover volume fraction for the cell model is found to be much smaller, 
rj w 0.2?y*. 

Working in the regime of saturated effective charges, we are left with in total only two independent input parameters: 
the volume fraction rj and the salt concentration na of the salt reservoir. Multiplying Z^f with As /a and the osmotic 
pressure with AttXbo 2 , as we have done in Eqs. ljA2|) . l|A5(l .(|4*|). we avoid As/a as an additional independent parameter. 
However, the colloid-colloid pair correlation function g{r) does not exhibit the same scaling behaviour as II m i cro so 
that it is important to precise which value of As /a has been used in the simulations. Unless otherwise specified, we 
have considered A^/a = 0.01, a reasonable value for colloidal systems in an aqueous environment (smaller values are 
also met in experiments). As for rj we have stopped our calculations at high volume fraction when the liquid started 
to solidify, while for na we have worked in between two extremes, na = 1.5 (high-salt regime) and na = 0.0 (no salt). 

Fig. compares the pressure prediction II m i cro of both models. The curves are based on Eq. (0J in combination 
with (|A3(I . (|A4|) . and (|A2(I . In the no-salt case both models have the same low-dilution behavior, but different limiting 
behavior in the presence of salt ions where we have an algebraic decay for the Jellium, but an exponential one for the 
PB cell model. The agreement between Jellium and PB-cell is excellent in the no-salt case, up to volume fractions 
around 0.1 which is remarkable in view of the differences in the effective charges. 

Before proceeding we wish to emphasize that many more effective-charge approaches can be found in literature, 
including DFT-based schemes (e.g. Monte Carlo solutions of the cell model (e.g. [UEl ): an d other mean-field 
models with various criteria for the effective charge (see reviews in [ill l28j | ) . We here have selected both the Jellium 
and the PB-cell models because they are complementary to each other and because they are probably the two most 
practical effective-charge models; they are sufficiently accurate, well-established in literature, quite simple and rather 
straightforward to implement (see appendix of |32lp. 




III. TEST I: MANY-BODY FORCES VERSUS YUKAWA FORCES 



We are now in the position to perform the first test. We recall that the interactions among colloids in a suspension 
are mediated by the microions and are therefore in origin complex many-body interactions. Integrating the microions 
out, the force on any given colloid depends on the positions of all other colloids in the system. It is not straightforward 
that this force can be written as a sum of pairs only. In principle, one has to sum over pairs, triple, ... and over 
all many-body configurations. To test how important the higher many-body contributions are, the true many-body 
forces have to be compared to the forces obtained by summing up the effective pair interactions |l4j . This is done in 
this section. 

We arbitrarily selected a few typical colloidal configurations from a MC simulation of 4000 Yukawa particles and 
solved the non-linear Poisson-Boltzmann equation in the region between the colloidal spheres of this colloidal con- 
figuration using our multi-centered Poisson-Boltzmann solver described in detail in |16|. As in [16l | we integrated 
the stress-tensor around each colloidal sphere to obtain the force acting on each colloid. These forces Fmb ~ which 
include all many-body forces acting on the individual colloid - can now be compared to the forces Fyu one obtains by 
summing pair- forces derived from Yukawa-potentials, Eq. Q, with effective PB-cell or Jellium parameters. We have 
evaluated the ratio of the magnitudes Fyu /Fmb and the (cosine) angles between the forces Fyu ■ Fmb/FyuFmb f° r 
each particle and have plotted the corresponding histograms, that have been averaged over a few configurations. The 
results of these calculations are shown in Fig. J2J) for the salt-free case and for the high salt case, na = 1.5. 

We see from Fig. (J2JA) that the ratio between the forces is very close to 1 and has a narrow delta-like distribution in 
the high salt case, as we should expect, since the Yukawa one-component picture has already proven to be applicable 
at such parameters. Interestingly, in the no salt case, the force ratio shows a broader distribution but is still quite 
peaked, although not centered at unity, see Fig. (J2JB). The forces based on the Yukawa one-component picture are 
therefore on average larger (cell) or smaller (Jellium) than the correct many-body forces. We may then anticipate 
that when used in the OCM, Jellium parameters will underestimate the structure, while cell data will lead to an 
overestimation (see Fig. 03). Despite these deficiencies, both models seem to reproduce the forces reasonably well, 
which is especially true for the angular distribution displayed in Fig. J2JC). An interesting result is the bump observed 
in Fig. J2JB) at low force ratios: here not any of the effective-charge models, but the Yukawa pair-force model as such 
fails dramatically, predicting a total force that is five times smaller in magnitude than the correct many-body force. 
We have investigated the local structure around the colloids feeling these forces. They are always located in regions of 
larger local density. Here, the local mean-distance between the colloidal particles is relatively short (compared to the 
suspension-wide mean distance) and the Yukawa pair potential is thus probed at rather short inter-particle distance. 
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FIG. 2: Histogram of ratio between the moduli of Yukawa pair forces Fyu and the correct many-body forces Fmb for A) high 
salt concentration (kci = 1.5) and for B) no added salt. A few particle configurations are obtained from MC simulations with 
4000 and 6000 particles at the volume fraction r\ — 0.0001. Fyu is evaluated with both effective parameters approximations, 
the PB cell and the Jellium. Fmb is obtained by solving the nonlinear Poisson-Boltzmann equation in the given colloidal 
configurations. The bare charge is such that ZXs/a = 50, which is large enough to effectively reach the saturation regime. 
We have also performed the calculations at larger volume fraction r\ = 0.01 and found very similar distributions as the ones 
shown here at r\ = 0.0001. C) The distribution of relative angles between the correct solution Fmb and the two effective pair 
pictures. All graphs are presented in a linear-log scale. 



However, at these distances expression Q is bound to fail, as we have indicated already in sec. 
this failure is now quantified in Fig. . 



The extent of 



IV. TEST II: THERMODYNAMIC CONSISTENCY AND THE KIRKWOOD BUFF-RELATION 



Having seen in the previous section that the distribution of forces felt in situ by the colloids in the mixture can in 
most cases be captured by a Yukawa effective pair potential, we now turn to the second test of these potentials and 
relate the compressibility of the suspension to its structure as computed within the OCM. The way to do so is to 
apply the Kirkwood-Buff relation introduced further below. 

n m i cro of both the PB-cell and the Jellium model is depicted in Fig. It provides an excellent approximation 
to the total pressure of salt-free suspensions: both models lead to a pressure that is in v ery good agreement with 
existing experimental data and primitive model simulations p3l Ej. see e.g. fA fill Eq. The Kirkwood-Buff 
relation (fTUfl now allows for a test of the effective potential chosen: u e g should lead to a colloid structure, embodied 
in a long- wavelength structure factor 5(0) computed within the OCM, that is compatible with the compressibility 
following from n m i cro . As the same effective parameters appear both in the Yukawa-forces leading to S(0) and in 
n m icro of eq. (0} this provides a critical consistency test of our effective-charge models. 

We will also report results at medium and high salt concentrations, where it does not seem possible to test the 
effective potential as severely as in the no salt case, but where it is still of interest to compare the different pressures 
introduced in this section. A similar idea has previously been worked out by Lobaskin et al. ^(| who checked the 
consistency of thermodynamic and structural properties of an asymmetric electrolyte containing macroions with 60 
elementary charges and monovalent counterions. Primitive-model results have been compared with different more 
approximate theories such as the cell model (within mean- field and even beyond). Finally, technical details may be 
found in Appendix 151 and IO 

We start this section with an outline of the theoretical background and present the actual results of our second test 
m section llVBl 
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A. Theoretical background 

1. The osmotic pressure of a charge-stabilized colloidal suspension 

Our problem is best treated in a semi-grand canonical ensemble. N c colloids are confined to a certain volume V 
which the small solvent molecules and additional micro-ions are free to leave; they are coupled to a reservoir fixing 
their chemical potentials. For the moment, we need not distinguish between the different types of small molecules, 
but may denote them collectively by an index s, while the colloidal degrees of freedom are referred to by the index c. 
The total internal energy of the system can then be written as U = U cc + U ss + U cs where the index combinations 
indicate over what degrees of freedom the bare interactions have to be summed. The force on particle i (i € c, s) then 
is Fi = — VjfJ and the total pressure of the colloidal suspension can be obtained using the virial 



P = pkT + — (yr i -F i ) (8) 




where the thermal average has to be taken over the colloidal as well as solvent /solute degrees of freedom. The first 
term, the ideal-gas term, is just the sum of densities of all components. 
Within the OCM, the pressure reads 



PckT+ — (Yri-Ff) (9) 




where the forces are obtained from the gradient of the effective potential and the summation runs over the N c colloids 
only. 

By definition, the effective potential u e e reproduces the same pair (colloid-colloid) distribution function g(r) as 
that of the original mixture, assuming a pair- wise Hamiltonian, see e.g. pj. However, there is no guarantee that 
the corresponding sum of effective pair forces coincides, for a given colloid i in a given colloidal configuration, with 
the micro-ion average of — Vjf7 which provides the true force. The comparison between both forces has been made 
in section HTT1 For the sake of the discussion, we assume here that they coincide. Eq. @ then uncovers one severe 
deficiency of the OCM picture : P ocm only provides a contribution to the total pressure given by Eq. (JSJ). In particular, 
in the limit of low salt the micro-ions are known to dominate the overall pressure of the suspension: that is the reason 
why Il m i cro of Fig. provides such an excellent approximation to the total pressure. This implies on the other 
hand that P OC m is a very poor approximation for P. The density derivatives of both pressures are however intimately 
connected, under all conditions, as we discuss in the following section (see eq. O). 

2. The Kirkwood-Buff relation 

The Henderson theorem [181 ] guarantees that every pair-correlation function g(r) can be uniquely associated with 
one pair-potential. As for the system considered here, this implies that for any given colloid density, there exists in 
principle one, and only one, effective pair-potential leading to a g(r) within the OCM that is in perfect agreement with 
the correct colloid-colloid pair correlation of the full multi-component system. Unfortunately, this does not guarantee 
that in such a multi-component system the higher-order correlation functions are equally well reproduced (see for 
example |19) where pair-, but not triplet correlations of a many-component colloidal systems could be reproduced 
within a simple pair-interaction picture). 

However, there is one important quantity one can still derive if the correct colloid-colloid g[r) is known: the osmotic 
isothermal compressibility x defined in eq. J3J). The Kirkwood-Buff relation relates this thermodynamic quantity to 
the infinite wavelength limit of the colloid-colloid structure factor S(q) [2(j 



dp c 

x - 



= 5(0) (10) 

T,salt 



d(3U 

where 5(0) is related to g(r) via 

S(0) = l + Pc I (g(r) - l)dr . (II 



Remarkably, Eq. (| 1 Of) is exact: nothing about the microion-microion correlations or the microion-colloid correla- 
tions needs to be known, but only the colloid-colloid pair correlations are required to compute the correct osmotic 
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compressibility of the full multi-component system. By contrast, the full compressibility of the system, defined as 
X" 1 = d(3P/d{p s + p c ) with P from Eq. @, depends on the pair correlations of all components Q. 

Connecting thermodynamic to structural information, Eq. Hl()|) is of central importance in the present work. In 
many colloidal suspensions the micro-ions determine the thermodynamics, while the macro-ionic degrees of freedom 
are more important for the structural properties of the system. In these cases Eq. <|10f) is also well suited to bridge 
the gap between the micro-ion and the macro- ion oriented viewpoints, thus providing a severe test for the quality of 
the effective potential. Indeed, though having a clear-cut definition, an effective potential is nevertheless a difficult 
object to compute and one often postulates its functional form, just as we have done in eq. (Q. With eq. (|10fl we can 
investigate a posteriori the relevance of the underlying ad-hoc assumptions made in deriving the effective potentials. 
Regarding the effective-charge concepts under scrutiny here, eq. (fTUf) describes a route how to connect u c s with 
n m icro- As the same effective parameters appear in both quantities, this relation also provides a consistency test of 
our effective-charge models. 



3. On the various contributions to the total pressure 



So far we have introduced three different pressures: P from which the osmotic pressure II follows, Poem, and finally 
-Pmicroj the latter quantity being for instance the pressure found in the cell or in the Jellium models, see eq. 
At this point, it seems worth discussing the non-trivial connection between these three pressure expressions, at the 
expense of introducing a fourth quantity: we define a "colloidal" pressure -P C oii from the difference between the total 
pressure P and P m i C ro 

P Pcoll ~t~ Pmicro- (1^) 

The pressure -P m i C ro can be considered as arising from so-called "volume terms" in the total free energy of the system, 
see e.g. HHE^. 

In the situation where the system is confined in a closed box of volume V (p c = N c /V), it can be shown that 



P = PckT+ i /^TrvffA + ^ <jf Pmicio (r)r-d 2 S^ , (13) 



The third term on the right hand side in (|13f) . for which the surface integral with normal oriented outward runs over 
the box confining the system, accounts for the direct coupling between colloids and micro-ions. It is precisely this 
quantity that the micro-ionic cell-model approaches aim at computing. In other words, this third contribution may 
be identified with P m i cr o in (|12(l so that remembering Eq. 0, we obtain P co i\ ~ Poem- However, a closed cell is not 
the most convenient configuration (in particular, the effective potential of interaction not only depends on the relative 
distance between two colloids, but also on the distance to the confining walls), and from a computational perspective, 
it is often more desirable to work with periodic boundary conditions systems. Unfortunately, Eq. ljT3)) which provides 
a physically transparent interpretation to P co ii m the confined case (close cell), breaks down with periodic boundary 
conditions. This failure is discussed at length in |25|, tog ether with the fact that in low salt (or no salt) conditions, 
one has P co n -C -Pmicro which implies P ~ P m icro (see [2(j for a related discussion). 

From the previous discussion, it appears that the connection between P ocm and P is not straightforward. The 
simplest relation between both quantities follows from the remark that the colloidal structure within the OCM is of 
course the same as in the original mixture (assuming the effective potential to be the correct one). Eas. l|10|) and i|ll|) 
then dictate that the compressibilities in both approaches coincide : 



dP 

dp c 



BP 

lJ1 ocm 

T,salt ®Pc 



(14) 

T,potcntial 



where it is crucial to compute the OCM compressibility at constant potential of interaction u e ff, i.e. discarding any 
density dependence of the effective potential. In a region of parameter space where the density dependence of u e g 
is absent or weak enough, the "salt" and "potential" subscripts in the partial derivatives of Eq. i|14[l correspond 
to the same constraint, since u e g then only depends on the salt chemical potential, besides relative distance. It is 
then possible to integrate Eq. (|14|l to obtain P ~ P OC m- Such a situation is met in the salt dominated regime to be 
discussed below. 



8 



1 TCOOOQOfa). 



X 



0.1 



0.01 



0.001 



0.01 



O S(0) 
X nrm (PB) 



X micro (PB) 

i/z;>b) 



Ka = 1.5 




0.0001 



0.001 



0.01 



0.0001 0.001 



11 



0.0001 



0.001 



0.01 



X 



T 



0.01 



0.01 



0.001 



Ka = 0.0 



r "q""o JO -O -O- . -Q — O. -QQ.QSi Q.O.Q OpB - ~~ 



0.001 . . oo 



_L 



0.0001 



0.001 



0.01 



FIG. 3: Compressibility vs volume fraction, for different salt conditions and applying different approximations as discussed in 
the text. Dashed line: Xmicro(PB), solid line: Xocm(PB) and empty circles: the colloid-colloid 5(0) from Fig. JjJJ. Dashed-dotted 
line shows 1/Zgg (PB), which is a valid approximation at all densities without salt, and holds at high densities only, when salt 
is present. 



B. Applying the Kirkwood-Buff relation 

We next apply eq. I|1(J|I to the case at hand and expect to find the following: In the no salt case where the total 
pressure is very close to P m i cr o we should obtain within the OCM a long-wavelength structure factor 5(0) compatible 
with dp c /d(pP m i cm ). Any difference between both quantities indicates a lack of consistency in the (somewhat un- 
controlled) procedure leading to the effective Yukawa potential JTJ chosen here. Such an inconsistency may take its 
roots in an incorrect computation of effective parameters (charge and screening length) or more fundamentally in the 
functional screened Coulomb form taken. 

We have calculated 5(0) within the OCM as a function of rj as described in detail in the appendix sectional (Monte- 
Carlo simulation) and ICl (integral-equation theory). Both methods lead to identical results, shown in Fig. © of the 
appendix. Such sets of 5(0) curves had to be computed twice, using once the PB cell model parameter and once the 
Jellium model parameter within the effective Yukawa pair-potential. Figures |3] (PB-cell) and 0] (Jellium) show the 
corresponding 5(0) curves. Each figure consists of three graphs, corresponding to Ka — 1.5, Ka — 0.5 and Ka = 0.0. 
In each graph, we compare 5(0) to the compressibility Xmicro following from the p c derivative of pressure P m i cr o (or 
equivalently II m i cro ). To distinguish between the two effective charge models used to compute II m i cro via eq. J1J, we 
introduce the notation Xmicro(PB) and Xmicro(Jell). Alternatively, one may compute the OCM pressure P ocm , defined 
in and obtains then the compressibility Xocm as the p c derivative of P cm- If PB-cell model (Jellium model) 
parameters have been used in the effective Yukawa potential, we denote the resulting compressibility by Xocm(PB) 
(Xocm(Jcll))- These quantities, and in addition the inverse of the saturated effective charges, are also graphed in 
Figure and E| 



1. No salt 



We observe that Xmicro, in the no salt case, is close to 5(0). The agreement seems to be slightly better for Jellium 
compared to cell parameters at low density, while the opposite holds at higher densities. This finding is consistent 
with the fact that the probability distribution shown in Fig. [5]B (where rj — 10~ 4 , a low value) is slightly more peaked 
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FIG. 4: Same as Fig. 10, but with effective parameters derived from the Jellium model. Dashed line: Xmicro(Jell), solid line: 
Xocm(Jell), empty circles: 5(0), dashed-dotted line: l/Zlg(Je\\) 



for the Jellium than for the cell. On the other hand, the inadequacy of Xocm for na = is expected : P ocm has here 
nothing to do with the total pressure P ; the resulting compressibility is about a factor of two larger than 5(0). This 
question is addressed in p5l | and for completeness, we adapt the argument in Appendix [DJ As can be seen in Figs. |3| 
and^l the no-salt compressibility is close to the inverse effective charge. This stems from the fact that at least within 
the Jellium, one has /3P m i cr o = Z c gp c for c s — and since the density dependence of Z c g is mild (logarithmic, see 
0), one has Xmicro — ■ Within the cell model, (iP m i cro — Z c gp c is no longer exact, but accurate enough to lead 

again to Xmicro — l/Z c g. 



2. Added salt 



In the first two figures in Figs|3]and0] (the ones with Ka > 0) we can see two regimes. At low volume fraction the 
OCM compressibility Xocm is a good approximation and the micro-ionic Xmicro fails, while at high enough volume 
fraction the reverse is true, Xmicro is good and Xocm fails. 

The derivative involved in the computation of Xocm includes the density dependence of the effective parameters. 
There are however situations where these parameters are virtually independent of p c (see Fig. |SJ. This is the case 
for example at low density or high enough salt concentration. Looking back at Eq. {TIJ, one realizes that in these 
cases the constant potential of interaction constraint coincides with that of constant c s . Then Eq. (|f 41) dictates that 
Xocm must be equal to 5(0). It is therefore not a surprise, nor a deep finding, that in Figs. [SJand^Jthe 5(0) coincides 
with Xocm in all those situations where the effective parameters have no p c dependence, i.e., at high enough salt 
concentration and at low density. 

Such an agreement notwithstanding allows us to compute the full pressure P when the effective potential is density 
independent. Remembering that at high salt content, the effective forces following from u e g are very close to their 
"exact" counterpart (see Fig. |2J), we expect P cm to provide there a good approximation to P, at least at not too 
high densities, i.e. in the salt-ion dominated screening regime when 77 <C 77* with rf from eq. 0. 

At high 77 the counter-ions dominate the screening. Thus, no matter how large Ka is, there is always a regime at 
high enough 77 where the salt-ions can be neglected and where the results become independent of Ka. That implies 
that curves differing in Ka must approach the same value at high enough r\. That can be observed, for example, in 
Fig. [SI A and Fig. of the appendix, where all curves show the same high 77 behavior. The same applies to Figures |3| 
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FIG. 5: A diagram showing the regions in parameter space of salt dominated screening and counter-ion dominated screening 
for both models. In the salt dominated screening regime (left from the dashed line) the OCM compressibility ^ cm is consistent 
with S(0), while in the counterion dominated screening regime (right from the dashed line for the cell model and right from 
the solid line for the Jellium) Xmicro is consistent with 5(0). On the solid line, rj is equal to rj* defined in eq. J7J, while on the 
dashed line r\ is equal to Q.2rf . 



and0] At high densities, all curves approach the zero-salt curves and go to 1/Z e g, which is the same for all na in this 
counter-ion dominated limit. On the other hand, for 77 — ?■ the compressibilities for ko ^ approach the ideal-gas 
limit. 

We have already remarked that for na = 0, Xmicro coincides with the OCM 5(0). Since, for large 77 (i.e. rj > 77*) all 
ko ^ curves must approach the na = curve, we expect to find that at high enough 77, Xmicro coincides with 5(0), 
since then the system is close to the salt- free regime |47| . This is indeed the case in Fig |3| (and to a lesser extent in 
Fig. 01 where only the trend is visible). This limiting behaviour seems to be a natural self-consistency requirement 
to impose to any effective potential. 

For the Jellium model we have already defined the crossover volume fraction 77* (eq. Q) separating the salt and 
counterion dominated screening. The dependence of rj* on the salt concentration na and on the volume fraction 77 is 
graphed in Fig. [S] (solid line). This graph is meant to summarize our findings of Fig. I^Jand^J For the Jellium model, 
the counterion dominated regime (77 > 77*) is to the right of the solid line. Here, the OCM 5(0) is found from Fig. 
0]to be consistent with Xmicro- The other regime (the salt-ion dominated regime) where we can observe from Fig. 0] 
that 5(0) is consistent with xocm is unfortunately not located to the left of the solid curve, but rather to the left of 
the dashed line empirically determined as 77* = O.277*. This means that in frame of the Jellium model there is a gap 
between the salt dominated region and counterion dominated region (77* < 77 < 77*). At these intermediate volume 
fractions neither Xmicro (Jell) nor Xocm(Jell) is consistent with 5(0). 

Such a gap is not found using parameters from the PB-cell model. For the PB cell model we have no analytical 
prediction for the value of the crossover volume fraction and have to read the values from the results in Fig. [3] 
Empirically we have found the crossover volume fraction of about 77 w O.277*, which is the same as the previously 
defined 77* for the Jellium model. This is plotted as a dashed line in Fig|5j For the PB-cell model, we may then 
summarize Fig. [3] as follows. For systems parameters lying on the right hand side of the dashed curve, 5(0) is 
consistent with Xmicro(PB), for those system parameters on the left hand side 5(0) is consistent with xocm(PB). 

We emphasize at this point that with salt at high 77, a better self-consistency is not necessarily synonymous with the 
fact that P m icro is in itself a better approximation for P. The situation is different for na = where we have the extra 
piece of knowledge that P ~ P m icro, that comes from comparison with experiments or primitive model computations 
(see also section W\ . 

After this discussion, we can summarize the results of our Kirkwood-Buff consistency check of the effective Yukawa 



11 




pair-potentials and the two effective charge models, as follows: 

• For na = 0, we find consistency essentially for both effective-charge models where the agreement is better for 
the Jellium model at low densities and for the PB-cell model at high densities. 

• For na = 0, the compressibility derived from P ocm is inadequate and not compatible with the OCM 5(0). The 
OCM approach with the Yukawa potential as effective pair-potential seems to produce correct results for the 
structure of the suspension, but not for the pressure. 

• The no-salt compressibility is close to the inverse effective charge. This offers a convenient way to estimate 
compressibilities Xmicro at na = 0. 

• For ko ^ and low volume fraction (77 < 77*), consistency is found between Xocm and 5(0), but not between 
Xmicro and 5(0). The OCM approach makes sense for calculating both the structure and the pressure. 

• For Ka ^ and high volume fraction (77 > rf for Jellium, 77 > 77* for PB-cell), we get back to the zero-salt case 
where 5(0) is consistent with Xmicro, but not with Xocm- 

• The value of r)* — 0.27/* has been empirically determined at the constant value of As/a = 0.01. This result 
might depend on the value of As/a, which has not been studied in scope of this paper. 

In cases where we have reason to assume that the OCM 5(0) is close to the correct x, we can of course estimate 
the pressure Ps(o) from integrating the OCM 5(0) 



^^ = 1 smn dp '- (15) 

To illustrate this idea, we have performed this integration in Fig. for na = and Ka = 1.5, using the PB-cell effective 
parameters. The resulting curves can then be compared to II m i cro from cq. 10} (again with the PB-cell parameters), 
which for Ka = we know to be an excellent approximation to the full pressure of the suspension. We have data 
for 5(0) starting at volume fraction 770 ~ 10~ 8 , which is used as the lower bound of integration in Eq. I] 15(1. A 
small constant is then added to the right hand side of Q15|). so that the result coincides with n m i cro at 770. The figure 
demonstrates that Ps(o) agrees well with n m i cro , surprisingly, not only for Ka — but also for Ka = 1.5. However, at 
very low 77, the agreement is excellent for Ka = 0, but not for Ka = 1.5, as one would expect from Fig. Q. Note that 
a hypothetical discrepancy between P m i C ro and Ps(o) here could not be considered as a lack of self consistency, since 
the contribution P co ii in P = P co ii + -Pmicro can be non-negligible. 
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FIG. 7: A) Osmotic pressure of a salt-free colloidal suspension according to the primitive model simulations of Ref. |43jj 
compared to the PB-cell and Jellium predictions of Eq. @. All data are computed for Zbarc A_e/a = 14.2 and As/a = 0.35. 
The effective charges are not saturated. If they were, one would obtain the PB cell pressure marked by the dashed-dotted 
line. B) Colloid-colloid pair-correlation functions from Linse's primitive model simulation (symbols) and from MC Yukawa 
simulations with effective PB-cell (solid lines) and Jellium (dashed lines) parameters, see also |4fjj |. From left to right, the 
packing fractions are r\ = 0.08, 0.04, 0.01, 0.0025, 0.00125. C) Compressibilities derived from the data in A) via Eq. © and 
from the colloid-colloid S(0) as indicated. 



V. TEST III: COMPARISON TO PRIMITIVE MODEL DATA 



P. Linse has carried out a full primitive model simulation of a salt-free charge-stabilized colloidal suspension and 
computed the pressure P as a function of volume fraction rj 0] , see also 0] . The highest colloidal charge considered 
in this work does not lead to fully saturated effective charges (^barc-Ve/a ~ 14.2). Fig. J7|A) shows the osmotic 
pressure as a function of r\ from |43( and compares it to P m icro (= n m i cro if na = 0). -P m icro is given by Eq. Q with 
either PB-cell or Jellium model data. For comparison, we also added the cell-model pressure one would obtain if the 
charges were saturated. This figure illustrates the quality of PB-cell and Jellium P m i cr o in the de-ionized limit, which 
has been repeatedly emphasized in the previous analysis. Fig. (JJJC) shows the corresponding compressibilities, and 
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as expected from Fig. (0A), it can be seen that Xmicro is very close to the correct primitive model compressibility. 
Figure QC) supports our findings of the previous sections: at low density, the Jellium effective potential leads to a 
5(0) that fares slightly better than its cell counterpart. The OCM cell and Jellium 5(0) respectively provide lower 
and upper bounds for the true compressibility, and correspondingly, respectively upper and lower bounds for the pair 
distribution function g(r) (see Fig. (0B)). This can be traced back to the force distribution seen in Fig. [5)3) where 
the Jellium effective potential leads to a slight underestimation of the "exact" force distribution, while the opposite 
holds for cell data. Somehow, considering an appropriate mean of Jellium and cell effective parameters would improve 
the quality of the predictions. A step forward in this direction has recently been made in |52(. 

VI. SUMMARY AND CONCLUSION 

The forces between colloids in a charge-stabilized colloidal suspension are commonly approximated by a sum of 
pair forces derived from the gradient of a screened Coulomb potential w c g with effective charge and screening length. 
A number of models can be found in literature on how to determine these effective parameters. We here computed 
them either from the Poisson-Boltzmann cell model supplemented with the Alexander et al recipe |3lT l35| , or from the 
renormalized Jellium model 0, H3 • The effective parameters are therefore not fitted but derived from a well defined 
although difficult to control procedure. It was our main motivation to assess the relevance of this procedure, be it in 
the Jellium or in the cell case. To this end, we have performed three tests. 

Test 1: We have performed in situ measures of effective forces in typical colloidal configurations, from the solution 
of a multi-center Poisson-Boltzmann theory (with N c — 4000 colloids). We are not aware of similar measures in the 
literature |l4j . The distribution of these forces has been compared to forces obtained from summing Yukawa effective 
forces. Under high salt conditions, the pair-wise and many body approaches gave very similar force distributions (the 
test was performed under conditions where cell and Jellium effective parameters are very close). In the no salt regime, 
the agreement is less quantitative but still quite good given that the system is there always strongly correlated, with 
strongly overlapping double- layers. Jellium(/cell) forces slightly underestimate(/overestimate) the "true" force. We 
also found situations where the Yukawa pair-force model fails dramatically, predicting a total force that is up to five 
times smaller in magnitude than the correct many-body force (bump observed in Fig. (|2lB) at low force ratios). 

Test 2: We have used the Kirkwood-Buff identity to check the consistency in the procedure leading to the effective 
Yukawa potential w e ff and the two effective charge models. On the whole, the consistency of u e s is remarkable given 
the simplicity of the underlying procedures. Our findings in detail have already been summarized at the end of 
section flVBI The most important results are: 

(i) In the salt-free case, we find consistency for both effective-charge models. The micro-ionic contributions Xmicro 
are consistent with the OCM 5(0). The Jellium potential is slightly better at low densities, but performs less than 
its cell counterpart at higher volume fraction. By contrast, the compressibility derived from P ocm is not compatible 
with the OCM 5(0) and consequently the OCM virial pressure P ocm is far from the true pressure P. 

(ii) With added salt, the P cm provides a good approximation to the total osmotic pressure if the volume-fraction 
T] is lower than the threshold value 0.2ry*, (77* of eq. Q). This is the case essentially if u e s is density independent, 
i.e., at low enough colloid density and/or high enough salt content. 

(iii) At Ka ^ and if 77 3> rj* , the system should recover a salt-free-like state: the cell rather than Jellium effective 
potential leads to a structure that is more consistent with the isothermal compressibility Xmicro- 

Test 3: We compared our data to primitive model (PM) data. PB-cell and Jellium both lead to a P m i C ro and 
Xmicro that is in excellent agreement with the corresponding PM pressure data. The OCM cell and Jellium 5(0) 
respectively provide lower and upper bounds for the true compressibility, and correspondingly, respectively upper and 
lower bounds for the pair distribution function g(r). 

Thinking in more practical terms, we finally recommend: 

(i) to use our empirical value 0.2r)* to find the regions where P cm can be taken as a good approximation of the true 
pressure and where thus the OCM picture is valid both with respect to the structure and to the thermodynamics, 

(ii) to always compute effective parameters applying both the Jellium and the PB model and then to take the 
results as an upper and lower bound. This latter remark is based not only on the results of the third test, but also on 
our observation that the correct forces J2JB) have always been located between the forces predicted by both models. 

(iii) to use the PB cell model to calculate the pressure -Pmicro(PB) at volume fractions above the threshold 77 > 0.2^*. 
This proves to be a remarkably consistent approximation (which does not necessarily mean correct !) of the true 
pressure P in this parameter region. 
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APPENDIX A: EFFECTIVE CHARGE MODELS 



1. The Jellium model 

The Jellium model 0, IS H(j assumes the effective charges of N c — 1 colloidal spheres to be smeared out in space to 
form a homogeneous background charge —pback = —ZbackPc- This background charge adds to the charge distributions 
of co- and counterions, p± = c s exp(=pe/30), in the radial PB equation, to be solved about one central colloidal 
particle. Deriving now Z e g by comparing the solution of the PB equation with the known far-field expression for 
the electrostatic potential <fi(r), one can again compute a colloidal background charge Zback which differs from the 
previous one. The whole procedure is iterated until self-consistency is achieved, that is, until Z c >&(Zback) — Zback- 
This consistency requirement was absent in the original formulation proposed in [2{j ; it is however an important 
ingredient to account for non- linear screening effects. Z e ff thus determined is used to calculate 

K 2 cS = 4ir\ B ^/z 2 s p 2 + 4c 2 (A1) 

which may be rewritten as 

4 4 ( Z c h\b \ 4 4 , . , 

n eS a = 13 — - — r) 1 + k a (A2) 

where 77 = 47rp c a 3 /3 is the colloid volume fraction while k 2 = 8ttXbc s is the squared inverse screening length in the 
reservoir. 



2. The Poisson-Boltzmann cell model 

The PB cell model rests on the observation that the repulsively interacting colloids arrange their positions 

such that each colloid has a region around it which is void from other colloids and which looks rather similar for 
different colloids. In other words, the Wigner-Seitz cells around two arbitrarily selected colloids are comparable in 
shape and volume. One now assumes that the total charge within each cell is exactly zero, that all cells have the 
same shape, and that one may approximate this shape such that it matches the symmetry of the colloid, i.e., spherical 
cells around spherical colloids. The cell radius R is chosen in consistency with the colloid volume fraction, and the 
PB equation within the cell is solved with appropriate boundary conditions at the cell edge and the colloid surface. 
Thus, through the finiteness of the cell plus the boundary conditions, the presence of all those colloids not inside the 
cell are taken into account. For a more detailed description of the cell model approximation see @- 

From the numerical solution of the PB equation, one obtains the electrostatic potential at the cell edge 4>r, and 
can now proceed to compute the effective screening parameter 

K 2 s a 2 = n 2 a 2 cosh<f> R (A3) 

if K 2 a 2 > and 

K 2 cS a 2 = ^ 2 exp(-(^) (A4) 

if K 2 a 2 — where p 2 appears in the PB equation V 2 = — p 2 e"^ of the salt-free case and is determined from the 
electro-neutrality condition. Then, the effective charge following Alexander and collaborator's recipe [3lJ is given by 

[mum 

Z cS ^- =7o/(w,rT 1/3 ) (A5) 
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FIG. 8: A) Effective saturated charges as a function of the volume fraction of a charge-stabilized colloidal suspension, for 
three different salinities, according to the Jellium model (symbols) and the PB-cell model (lines). We can discriminate the salt 
dominated regime r\ < rf and the counterion dominated regime r\ > 77* . B) The radial distribution function plotted versus 
the scaled distance d/d nn for no-salt systems at four different volume fractions: 0.00125 (the lowest peaks), 0.005, 0.02 and 
0.08 (the highest peaks). There are some differences, but in the first approximation the curves superimpose. Inset: The ratio 
l/(K e ffdnn) of the double layer thickness (derived from the cell approximation) and the mean distance between particles as a 
function of the colloid volume fraction r\ at three different salt concentrations. 



where 70 = taring in the salt case and 70 = 1 in the no-salt case, and where the function f(x,y) is given by 

fix, y) = — \(x 2 y — 1) sinhfrcy — x) + x(y — 1) cosh(a;?/ — x)\ . (A6) 
x 

A simple approximation valid for large bare charges is given in |35| . We emphasize here that following such a route 
to define Z e g and K e ff, a "natural" relation such as that embodied in Eq. (|A1I) is lost, except in certain particular 
limits (low density, or high density, or low charge). We also note that it has been shown recently that the cell 
model effective charge is accurately reproduced by a dynamical rule which defines the condensed microions through 
a bound on their total energy [3r3 |. a criterion that may also be considered when mean-field breaks down. 

Fig. |JSJ serves to discuss and compare the two effective charge models in the (na, ^-parameter space. Fig. |[8JA) 
compares the volume-fraction dependence of the saturated effective charges. The change from salt-ion- to counter- 
ion-dominated screening (occurring at 77 ~ 77* ) can be recognized from the onset of a ^-dependence of the effective 
charges. For the no-salt case screening is always due to the colloidal counter-ions, and, indeed, both models predict a 
strong ^-dependence of the effective charges even in the limit rj — > 0. At low 77, a suspension of charged colloids with 
no extra salt will always be correlated regardless of how low a volume fraction is considered. The reason is that the 
thickness of the double layer then grows faster than the mean-distance between the particles, that is, the ratio of the 

— 1/3 

double-layer thickness 1 / K e ff and the mean distance d nn = Pc grows with decreasing volume fraction, as opposed 
to the case with external salt where this ratio decreases (see and the inset of Fig. |8jB). 

Fig. @AJ_demonstrates that the two effective charge models agree inasmuch as the 77 — > limiting behavior is 
concerned |39j, but disagree in the opposite limit with the Jellium model generally predicting an earlier change from 
the salt-ion- to the counterion-dominated screening and smaller effective charges. For high 77, the curves for all three 
values of na must ultimately converge for each model when the contribution of salt-ions to the screening ceases to 
be significant. A rather special feature of the Jellium model is the pronounced minimum in Z^S at intermediate 
volume fractions observed for all values of kcl, something that is not present in the PB cell model for the salinities 
investigated, but that would be observed at lower salt. 



APPENDIX B: NUMERICAL PROCEDURES 



Technical details for the calculation of the effective charges and screening parameters are given in |7], |30( . Within 
the OCM picture, we have performed Monte-Carlo simulations, typically with 10000 particles in a simulation box of 
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FIG. 9: Comparison of the large wavelength limit of the structure factor 5(0) as calculated by Monte Carlo simulations and by 
OZ integral equations method. The cell model approximation is used to obtain the effective parameters. The symbols represent 
the MC results, the solid lines OZ values. 



side length Lf, ox applying periodic boundary conditions. We have carried out 5 • 10 6 MC cycles to reach thermal 
equilibrium and another 2 • 10 7 cycles for the averages. The pressure has been obtained from Eq. ©, the structure 
factor directly from 

= ^ ( (E C0S (S^)) + (j2 sin(q-r^ ^ (Bl) 

and 5(0) has then been approximated by S(q m i n ) where q m in — ^n/Lbox |48| . Special care has been taken that the 
value converges with the system size. As an independent check we computed the static structure factor S(q) from the 
g(r) by Fast Fourier Transform and checked that the value at q m i n was the same. 

Much faster than MC simulations are structure calculations using the Ornstein-Zernike (OZ) equation |4pj . see 
appendix |0 This integral equation has been solved using the well-tested Rogers- Young closure j41j. Computing thus 
g(r), we obtain 5(0) from the integration in Eq. I|ll|) supplemented with finite size scaling, and the pressure from 

pPoom = Pc-§J rg(r)(3u' cS (r)dr . (B2) 

To assess the validity of both Monte Carlo and integral equation routes, we show in Fig. © the results of our 5(0) 
calculations in which the PB-cell effective parameters have been taken in the Yukawa potential. The results of the 
MC simulations compare favorably with the solutions of the OZ equation for all values of na considered. Equally 
good agreement between the results of both methods was found for all other curves presented in this work, but in 
order not to overload the graphs, we have shown only the OZ results. 



APPENDIX C: INTEGRAL EQUATIONS THEORY 

The structure of a fluid can be obtained by experimental techniques, computer simulations, or by solving numerically 
the Ornstein-Zernike (OZ) equation. The inhomogeneous OZ equation for a mono-component fluid is given by |40| 

h(r 1 ,r 2 ) = c(ri,r 2 ) + / dr 3 c (ri, r 3 ) p (r 3 ) h (r x , r 3 ) , (CI) 
Jv 

where h (ri, r2) and c (ri, r 2 ) are the total and direct correlation functions, respectively, between a particle located at 
ri and a particle located at r 2 . p (r 3 ) is the local density of particles in the system. In an homogeneous and isotropic 
system the total and direct correlation functions depend only on the relative distance between particles and the local 
density takes the average value p(r 3 ) = p, where p is the mean density of particles. Then, equation (|C1|) reduces to 



h{r)=c{r)+p [ dv'c{\v - r'Q h{\v - r'|) 
Jv 



(C2) 
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The total correlation function is related to the local structure of the system by means of the relation h (r) = g (r) — 1, 
where g (r) is the radial distribution function. Also, h(r) is connected with the structure factor, S(q), through the 
relation S (g) = 1 + h (q) , where h (q) is the Fourier transform of h (r) . The structure factor can be measured, for 
example, by light scattering experiments. 

Eq. i|C2|) is an integral equation with two functions, h (r) and c (r), which needs an additional relation between the 
total and direct correlation functions to close the set of equations. The general closure relation for the equation (IC2I) 
is given by 01 



h (r) = exp [-(3u (r) + h (r) - c(r) + B (r)] - 1, 



(C3) 



where f3u (r) is the pair potential between particles and the function B (r) is the so-called bridge function, which 
depends on particle density and, in general, is unknown. There are many approximations to the b ridg e function, 
such as Percus-Yevick (PY), hyper-netted chain (HNC) and Rogers- Young (RY) closure relations [40l l4l| It is well- 
known that the HNC and RY relations work well when the interaction between particles is only repulsive beyond the 
hard-core interaction. The RY approximation |4l| is given by 



h (r) = exp [— (3u{r)\ 



1 



exp[(fe(r)-c(r))/(r)]-l 

m 



-i, 



(C4) 



where the function / (r) = 1 — exp (—or). The RY closure is a mixture between the PY and HNC closure relations. 
For example, when a = Eq. (|C4(I reduces to PY approximation. As a increases, / (r) approaches 1, and Eq. (|C4|) 
reduces to HNC approximation. The mixture parameter, a, is fixed by demanding that the isothermal compressibility 
from the virial route, 1 — pc(q = 0) = S(q = 0), and from the compressibility route (see Eq. ©), are both equal. 
In the case of density-dependent potential of interaction, as is the case here, it is important to compute the density 
derivatives leading to the above compressibilities at constant potential of interaction. 

For a given density p the integral equation (|C2jl . together with the RY closure relation, is numerically solved by 
converting the OZ equation in an algebraic equation. This is done by taking the Fourier transform of Eq. I|C2|) . 



h(q) 



c(q) 



1 - pc{q) ' 



(C5) 



Then, by applying the inverse Fourier transform to Eq. I|C5(I we get the desired solution. However, for strongly 
interacting systems, as in our case, the direct application of Eq. l)C5j) gives noisy solutions. Then, we divide the 
density p in small steps of size Ap. For each given sub-density we solve iteratively the OZ by using a five-parameter 
version of the Ng-method 0] until the desired density is reached. At each step of the iteration, the pair distribution 
function is determined and used together with the RY closure relation in order to compute the new direct correlation 
function, c(r). To ensure rapid convergence, the value of c(r) at the previous step is taken as an initial guess. 



APPENDIX D: A SIMPLE APPROXIMATION FOR THE OCM PRESSURE 

We concentrate on the salt- free case. With a potential given by Eq. a straightforward computation of P OC m as 
given in (JBU leads to (lll2l| 

PP~» - -'f r 9{r) d -^r^r (Dl) 

6 Jr=2a dr 



^Pl z ls^B \, , (KcSa) 2 \ ( pf 



^cff 



3(1 + K c ff a) 2 J 6 



[g(r) - 1] (1 + KcS r) (3u eS (r) d 3 r 



where the ideal gas contribution present in l|B2fl is a small quantity (whenever Z e g 3> 1 which is the case for highly 
charged colloids where Z e g is typically of order 10a/ Xb) and has been discarded. In IjDlJl . the dominant term is the 
first one, arising from the long-range behavior of the pair correlation function (g — > 1 at large distances) . In this term, 
the curly brackets may be safely approximated by 1 since at low densities, K c ga <^ 1. Therefore 

PP ocm - ^#t^. (D2) 

Within the jellium, we have K 2 ff = AnX^Z^pc and /3P m i CTO = Z e gp c (that are approximately correct within the cell 
model), so that 

ftPocm — -^Z cS p c = -/3P mi cro- (D3) 
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Hence the factor of 2 between Xocm and x m i cr o, roughly observed in Figs. [3] and 0] in the na = case. The fact 
that Xocm 7^ 5*(0), or in other words, that the virial and the compressibility routes do not coincide, is a well known 
deficiency of density-dependent pair potentials, see e.g. f° r a general discussion. Remaining at the OCM level, 
improving the reliability of the virial route - which takes P ocm to approximate the total pressure - is possible by 
formally including a density derivative of the effective potential into the forces from which the virial is computed, see 
e.g. sections 4.2 of both references pj or j^jj. We do not follow such a route here, since the purpose is not to test 
improvements of the virial route by changing P ocm while neglecting P m i cr o, but to test the internal consistency of a 
procedure (cell or jellium) that leads both to P m i C ro and to the effective potential u e g. 

Equation i|D2|) also offers a clear illustration of Kirkwood-Buff identity i|14|) . To compute the right hand-side of 
(|14f) . we have to fix Z g and n e g while computing the derivative, so that 



dP oc 



dp c 



47r p c Z 2 ff X B 



T. potential 



(D4) 



On the other hand, since P ~ P n 



dP 

dp c 



dZ cfi 



T.salt 



dp c 



(D5) 



where the second term on the right hand side is negligible, as already argued in section llVl We consequently see that 
Eq. (|14|) is fulfilled, within the approximations invoked. 
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